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Abstract. Acyclic and cyclic orientations of an undirected graph have 
been widely studied for their importance: an orientation is acyclic if it 
assigns a direction to each edge so as to obtain a directed acyclic graph 
(DAG) with the same vertex set; it is cyclic otherwise. As far as we know, 
only the enumeration of acyclic orientations has been addressed in the 
literature. In this paper, we pose the problem of efficiently enumerating 
all the cyclic orientations of an undirected connected graph with n ver¬ 
tices and m edges, observing that it cannot be solved using algorithmic 
techniques previously employed for enumerating acyclic orientations. We 
show that the problem is of independent interest from both combinato¬ 
rial and algorithmic points of view, and that each cyclic orientation can 
be listed with 0(m) delay time. Space usage is 0(m) with an additional 
setup cost of O(n^) time before the enumeration begins, or 0{mn) with 
a setup cost of 0(m) time. 

1 Introduction 

Given an undirected graph G{V,E) with n = \V\ vertices and m = \E\ edges, 
an orientation transforms G into a directed graph G by assigniM a dir ection 
to each edge. That is, an orientation of G is the directed graph G {V,Ej) such 
that the vertex set V is the same as G, and the edge set E is an orientation 
of E: exactly one direction between ( ' and (u, u) € holds for any 



undirected edge {u, u} €An orientation G is acyclic when G does not contain 
any directed cycles, ^ G is a directed acyclic graph (DAG); otherwise we say 
that the orientation G is cyclic. 

Acyclic orientations of undirected graphs have been studied in depth. Many 
results concern the number of such orientations: Stanley m shows how, given 
a graph G and its chromatic number x, the number of acyclic orientations of G 
can be computed by using the chromatic polynomial (a special case of Tutte’s 
polynomial). Other results rely on the so called acyclic orientation game: Alon 
et al. |T] inquire about the number of edges that have to be examined in order 
to identify an acyclic orientation of a random graph G; Pikhurko m gives an 
upper bound on this number of edges in general graphs. The counting problem 
is known to be ^P-complete [8] and enumeration algorithms that list all the 
acyclic orientations of a graph are given in [2] and m- 

We consider cyclic orientations, which have been also studied from many 
points of view. Counting them is co-#P-complete [8]. In Fisher et al. [5], given 


a graph G and an acyclic orientation of it, the number of dependent edges, i.e. 
edges generating a cycle if reversed, has been studied. This number of edges 
implicitly gives a hint on the number of cyclic orientations in a graph. In m 
an algorithm has been given to make inference about causal structure in cyclic 
graphs. In m directed cyclic graphs are used to model economic processes, and 
an algorithm is given to characterize conditional independence constraints of 
these processes. 

In this paper we address the problem of enumerating all the cyclic orienta¬ 
tions of a graph. 

Problem 1. Enumerating the set of all the directed graphs ^ that are cyclic 
orientations of an undirected graph G. 

We analyze the cost of an enumeration algorithm for Problem [T] in terms of 
its setup cost, meant as the initialization time before the algorithm is able to lists 
the solutions, and its delay cost, which is the worst-case time between any two 
consecutively enumerated solutions (e.g.^H). We are interested in algorithms 
with guaranteed 0{m) delay, where the O notation ignores polylogs. 

A naive solution to Problem [T] uses the fact that enumeration algorithms 
exist for listing acyclic orientations mM- It enumerates the cyclic orientations 
by difference, namely, by enumerating all the 2™ possible edge orientations and 
removing the a acyclic ones. However, this solution does not guarantee any 
polynomial delay, as the number /? = 2™ — a if cyclic orientations can be much 
larger or much smaller than the number a of acyclic ones. For example, a tree 
with m edges has a = 2™ and /3 = 0. On the other extreme of the situation, 
we have cliques. An oriented clique is also called a tournament, and a transitive 
tournament is a tournament with no cycles. A clique of n nodes (and m = 
edges) can generate 2™ different tournaments, out of which exactly a = n! will 
be transitive tournaments [5]. As 2™ grows faster than n!, we have that the ratio 
a/P tends to 0 for increasing n, where /3 = 2"* — a. 

To the best of our knowledge, an enumeration algorithm for Problem [1] with 
guaranteed 0{m) delay is still missing. We provide such an algorithm in this 
paper, namely, listing each cyclic orientation with 0(rn) delay time. Space usage 
is 0{m) memory cells with a setup cost of O(n^) time, or 0{mn) memory cells 
with a setup cost of 0{m) time. Interestingly, Problem [T] reveals to be a rich 
source for enumerations techniques, and our solution offers new combinatorial 
and algorithmic techniques when compared to previous work on the enumeration 
of acyclic orientations las]. 

In the following, for the sake of clarity, we will call edges the elements of 
E (undirected graph) and arcs the elements of 1^ (directed graph). We will 
assume that the graph in input G is connected and we will denote as n = \V\ 
and m = \E\ respectively its number of nodes and edges. 

The paper is organized as follows. Section [2] gives an overview of our enumer¬ 
ation algorithm. Section [3] describes the initialization steps, and shows how to 
reduce the problem from the input graph to a suitable multi-graph that guaran¬ 
tees to have a chordless cycle (hole) of logarithmic size. The latter is crucial to 
obtain the claimed delay. Section |4] shows how to enumerate in the multigraph 
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and obtain the cyclic orientations for the input graph. Section describes how 
to absorb the setup cost using more space. Finally, some conclusions are drawn 
in Section [HI 


2 Algorithm overview 


The intuition behind our algorithm for an undirected connected graph G(Vg , Eq) 
is the following one. Suppose that G is cyclic, otherwise there are no cyclic 
orientations. Consider one cj^l^G {Vc ^ Ec) in G: we can orient its edges in two 
ways so that the resulting G is a directed cycle. At this point, any orientation 
of the remaining edges, e.g. those in Eq \ Ec, will give a cyclic orientation of G. 
Thus, the interesting cases are when the resulting is not a directed cycle. 

The idea is first to generate all possible orientations of the edges in Eg \ Eg , 
and then assign some suitable orientations to the edges in Eg- This guarantees 
that we have at least two solutions for_^ach orientation of Eg \ Eg, namely, 
setting the orientation of Eg so that G is one of the two possible direjded 
cycles. Yet this is not enough as we could have a cyclic orientation even if G is 
not a directed cycle. 

Therefore we must consider some cases. One easy case is that the orientation 
of Eg \ Eg already produces a directed cycle: any orientation of Ec will give a 
cyclic orientation of G. Another easy case, as seen above, is for the two orien¬ 
tations of Eg such that G is a directed cycle: any orientation of Eg \ Eg will 
give a cyclic orientation of G. It remains the case when the orientations of both 
Eg \ Eg and Eg are individually acyclic: when put together, we might have or 
not a directed cycle in the resulting orientation of G. To deal with the latter 
case, we need to “massage” G and transform it into a multigraph as follows. We 
refer the reader to Algorithm [TJ 

Algorithm setup is performed as described in Section [3] We preprocess G 
with dead-ends removal and edge chain compression. The result is an undirected 
connected multigraph M{Vm, Em), where the edges are labeled as simple and 
chain. After that we find a chordless cycle of logarithmic size G in M, and 
remove Ec from Em, obtaining the labeled multigraph M'{Vm,E)^), where 
E'm = Em \ Ec- 

Enumerating cyclic orientations described in Section |4] exploits the prop¬ 
erty (which we will show later) that finding cyclic orientations of G corresponds 
to hnding particular orientations in M', called extended cyclic orientations, and 
of G, called legal orientations. In the for loop, these orientations of M' and G 
are enumerated so as to find all the cyclic orientations of G. As we will see for 
the latter task, it is important to have G of logarithmic size to guarantee our 
claimed delay. 


^ This will actually be a chordless cycle of logarithmic size (called log-hole). 
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Algorithm 1: Returning all the cyclic orientations of G 
Input: An undirected connected graph G(V, E) 

Output: All the cyclic orientations u {V, E) 

Algorithm setup (Section 151): 

Remove dead-ends (nodes of degree 1) recursively from G 

M{Vm, Em) t— replace G’s maximal paths of degree-2 nodes with chain edges 

G{Vc, Ec) t- a log-hole of M 

M'(Vm, E'j^) •<— delete the edges of C from M, i.e. E';^ = Em \ Ec 

Enumerate cyclic orientations (Section HI): 

for each extended orientation A?' of multigraph M' do 

for rach leg al o rientation of log-hole: G (see Algorithm ID) dp 
M"{Vm,E") t— combine M' and C, where E" = E'm'JEc 
Output each of the cyclic orientations C? of G corresponding to A?" 


3 Algorithm Setup 

3.1 Reducing the problem to extended cyclic orientations 

In the following we show how to reduce Problem [T] to an extended version that 
allows us to neglect dead-ends and chains. 

Dead-end removal. Given an undirected graph G(R, E), a dead-end is a node 
of degree 1. We recursively remove all dead-ends^ so that all the surviving nodes 
have degree 2 or greater. By removing these nodes and computing the cyclic 
orientations in the cleaned graph, we can still generate solutions for the original 
graph by using both orientations of the unique incident edge to each dead-end, 
as emphasized by the following lemma, whose proof is straightforward. 

Lemma 1. Let G be a graph, u a dead-end and {u,v} its unique incident edge. 
Let G' be the graph G without u and the edge {u, The set of all the cyelic 
orientations of Gis composed by the orientations G'{V' U |u|, 1^' 
and G'{V'U{u},E'U{{v,u)}), for all the cyclic orientations G'(V',E') of G'. 


The dead-end removal can be done by performing a DFS recursive traversal of 
G, starting from an arbitrary node x. Every time a node of degree 1 is visited, 
it is removed from the graph. When the recursion ends in a node, the latter is 
removed if all of its neighbors have been removed except one (which is its parent 
in the DFS tree). Finally, when the traversal ends, it might be that the node 
from which we started has degree 1. To complete the process, if x has now degree 
1, we remove it from the graph. The DFS of G and the removal of its dead-ends 
can be done in 0{m). 

The rationale for removing dead-ends is to have shorter cycles: for example, 
consider a “necklace” graph with n nodes, for n even, such that n/2 nodes form a 
cycle, and the remaining n/2 nodes have degree 1 and are attached to one of the 
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nodes in the cycle, such that each node in the cycle has degree 3 and is connected 
to one node of degree 1. With the removal of dead-ends, the cycle has only nodes 
of degree 2 and can be compressed as discussed in the next paragraph. 

Chain compression. It consists in finding all the maximal paths vi,... ,Vk 
where Vi has degree 2 (with 2 < i < k — 1), and replacing each of them, called 
chain, with just one edge, called chain edge. It is easy to see that this task 
can be accomplished in 0{m) time by traversing the graph G in a DFS fashion 
from a node of degree > 3. The output is an undirected connected multigraph 
M{Vm,Em), where Vm C V are the nodes of V whose degree is > 3, and Em 
are the chain edges plus all the edges in E which are not part of a chain@ The 
latter ones are called simple edges to distinguish them from the chain edges. In 
the rest of the paper, M will be seen as a multigraph where \Vm\ > 4 and each 
of the edges has a label in {simple, chain}, since it might contain parallel edges 
or loops. For this, we define the concept of extended orientation as follows. 

Definition 1 (Extended Orientation). For a multigraph M{Vm, Em) hav¬ 
ing self lop ps and edges lab e lec) as simple and chain, an extended orientation 

{Vm,Em) is an orientationof its edge set Em'- for any simple edge {u, ?;}, 
exactly one direction between {u,v) and {v,u) holds; for any chain edge {u,u}, 
either is broken, ^or exactly one direction between {u,v) and {v,u) holds. A di¬ 
rected cycle in M cannot contain a broken edge. 

Broken edges correspond to chain edges that, when expanded as edges of G, 
do not have an orientation as a directed path. This means that they cannot be 
traversed in either direction. Note that this situation cannot happen for simple 
edges. The following lemma holds. 

Lemma 2. If we have an algorithm that list all the extended cyclic orientations 
of M(Vm, Em) with delay f{\EM\), for some / : IR —>■ IR, then we have an 
algorithm that lists all the acyclic orientations of the graph G{V, E) with delay 
0{f{\EM\) + \E\). 

Proof. For each extended cyclic orientations A? we return a set S of cyclic orien¬ 
tations of G: any edge e of A? maintain^he same direction specified by A? in all 
the solutions in S; for each chain c oi M, we consider the edges corresponding 
to c in G, say ei, 62 ,..., Ch'- if c has a direction in A?, the same direction of c 
is assigned to all the edges Cj in all the solutions in 5; if c has no direction 
assigned, i.e. broken, we have to consider all the possible 2^—2 ways of making 
the path ei, 62 ,..., e/j broken (these are all the possible ways of directing the 
edges except the only two directing a path). All the solutions in S differ for the 
way they replace the chain edges. 

Getting extended cyclic orientations in /(|Am|) delay, iterating over all the 
chain edges c, and iterating over all the corresponding edges of c assigning the 
specified directions as explained above, we return acyclic orientations of the 
graph G(y,E) with delay 0(/(|i?M|) + lAl). □ 

The degenerate case of M with < 3 nodes can be handled separately. 
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Lemma [5] allows us to concentrate on extended cyclic orientations of the 
labeled multigraph M rather than on cyclic orientations of G. Conceptually, we 
have to assign binary values (the orientation) to simple edges and ternary values 
(the orientation or broken) to chain edges. If we complicate the problem on 
one side by introducing these multigraphs with chain edges, we have a relevant 
benefit on the other side, as shown next. 


3.2 Logarithmically bounded hole 

Given the labeled multigraph obtained in Section IXTl namely M{Vm, Em), we 
perform the following two steps. 

1. Finding a log-hole. Find a logarithmically bounded hole (hereafter, log- 
hole) C{Vc, Ec) in M{Vm,Em)' it is a chordless cycle whose length is either 
the mrth of the graph (i.e. the length of its shortest cycle) or this length plus 
one@ 

2. Removing the log-hole. Remove the edges in Ec from M, obtaining 
M'{Vm,E'j^), where E'j^ = Em - Ec- 

Properties of the log-hole. Since M is a multigraph with self-loops, a log- 
hole C{Vc, Ec) of M can potentially be a self-loop. In any case, the following 
well-known result holds. 

Lemma 3 (Logarithmic girth | [3]|4] 1. Let G{V,E) be a graph in which every 
node has degree at least 3. The girth of G is < 2[log|R|]. 

Lemma [3] means that the log-hole G oi M has length at most 2 [log \ Vm 11+1, 
thus motivating our terminology. 

The log-hole G can be found by finding the girth, that is performing a BFS 
on each node of the multigraph M to identify the shortest cycle that contains 
that node, in time 0 (|Vm| • \Em\)- By applying the algorithm in [6], which easily 
extends to multigraphs, we compute G in time 0{\Vm\’^)'- in this case, if chords 
are present in the found G, in time 0(log \ Vm\) we can check whether G includes 
a smaller cycle and redefine G accordingly. 

4 Enumerating cyclic orientations 

We now want to list all the cyclic orientations of the input graph G. By LemmaH] 
this is equivalent to listing the extended cyclic orientations of the correspond¬ 
ing labeled multigraph M{Vm, Em) obtained from G by dead-end removal and 
chain compression. We now show that the latter task can be done by suitably 
combining some orientations from the labeled multigraph M'{Vm, E'j^) and the 
log-hole G{Vc, Ec) using an algorithm that is organized as follows. 

® Minimum cycle means the cycle having minimum number of edges (i.e. a self loop), 
where chain edges count just one like normal edges. 
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1. Finding extended prientations. Enumerate all extended orientations (not 

necessarily cyclic) M' of the multigraph M'. _^ 

2. Putting back the log^hole.For each listedM'(VM 5 -E'^M), consider all_the 
extended orientations C {Vc,l^c) of the log-hole C such that E'm ^ Ec 
contains a directed cycle, and obtain the extended cyclic orientations for the 
multigraph M. 

Finding extended orientations. This is now an easy task. For each edge 
{u,ri} in E'j^ that is labeled as simple, both the directions (u,v) and {v,u) 
can be assigned; if {u,z;} is labeled as chain, the directions (u,v) and (v,u), 
and broken can be assigned. Each combination of these decisions produces an 
extended orientation of M'{Vm,Em)- If there are s simple edges and b broken 
edges in M', where s -I- to = this generates all possible 2®3^ extended 

orientations. Each of them can be easily listed in 0{\E'j^^\) delay (actually less, 
but this is not the dominant cost). 


Putting back the log-hole. For each listed M' we have to decide how to put 
back the edges of the cycle C, namely, how to find the orientations of C that 
create directed cycles. 

Definition 2. Given the cycle C{Vc, Ec) and ivf'(V m, ^m)) we call legal ori¬ 
entation cf (Vn. ^n) anv extended orientation ofC such that the resulting multi- 



The two following cases are possible. 

1. A?' is cyclic. In this case each edge in Ec can receive any direction, including 
broken if the edge is a chain edge: each combination of these assignments 
wi ll produce a legal orientation that will be output. 

2. M' is acyclic. Since C is a cycle, there are at least two legal orientations 
obtained by orienting C as a directed cycle clockwise and counterclockwise. 
Moreover, ^ding just an oriented subset of edges D C C to can create 
a cycle in M': in this case, any orientation of the remaining edges oi C\D 
(including broken for chain edges) will clearly produce a legal orientation. 

While the first case is immediate, the second case has to efficiently deal with 
the following problem. 

I^oblern.2. Given M' acyclic and cycle C, enumerate all the legal orientations 


d{Vc,Ec) of C. 


In order to solve Problem [21 we exploit the properties of C. In particular, 
we compute the reachability matrix R among all the nodes in Vc, that is, for 
each pair rt, v of nodes in Vc, R{u^ v) is 1 if it can reach v in the starting M, 
0 otherwise. We say that R is eyelie whether there exists a pair i, j such that 
R{hj) = Rijii) = 1- This step can be done by performing a BFS in M' from 
each node in Vc- by ?? we have iVcl < 2|’log|VMn + 1; so the cost is 
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0{\E'j^\ ■ log 114^1) time. Deciding the orientation of the edges and the chain 
edges in Ec is done with a ternary partition of the search space. Namely, for 
each edge {u,r’} in Ec, if {m,^} is simple we try the two possible direction 
assignments, while if it is a chain we also try the broken assignment. In order to 
avoid dead-end recursive calls, after each assignment we update the reachability 
matrix R and perform the recursive call only if this partial direction assignment 
will produce at least a solution: both the update of R and the dead-end check 
can be done in 0(log^ I^mI) (that is, the size of R). 

Scheme for legal orientations. The steps are shown by Algorithmic] At the 
beginning the reachability matrix R is computed and is passed to the recursive 
routine LegalOrientations. At each step, C ' is the partial legal orientation to 
be completed and / is the set of broken edges declared so far. Also, j is the index 
of the next edge {cj, c^+i} of the cycle C, with 1 < j < h (we assume Ch+i = Ci 
to close the cycle): if j = h -|- 1 then all the edges of C have been considered 
and we output the solution u' together with the list / of broken edges in cf'. 
Each time the procedure is called we guarantee that the reachability matrix R 
is updated. 

Let {u, u} be the next edge of C to be considered: for each possible direction 
assignment (u, v) or (u, u) of this edge, we have to decide whether we will be able 
to complete the solution considering this assignment. This is done by trying to 
add the arc to the current solution. If there is already a cycle, clearly we can com¬ 
plete the solution. Otherwise, we perform a reachability check on {cj+i,... Ch+i}'. 
it is still possible to create a directed cycle if and only if any two of the nodes 
in {cj+i,... Cfi+i}, say c/ and Cg satisfy R{cf,Cg) = 1 or R{cg,Cf) = 1. This 
condition guarantees that a cycle will be created in the next calls, since we know 
there are edges in C between c/ and Cg that can be oriented suitably. Finally, 
when {u, u} is a chain, the broken assignment is also considered: R does npt 
need to be updated as the broken edge does not change the reachability oiM'. 

The reachability and cyclicity checks are done by updating and checking the 
reachability matrix R (and restoring R when needed). Updating R when adding 
an arc (u, v) corresponds to making v, and all nodes reachable from v, reachable 
from u and nodes that can reach u. This can be done by simply performing an or 
between the corresponding rows in time 0(log^ |Um|)j since i? is \C\ x \C\. The 
reachability check can be done in 0([o^ \Vm\) time. The cyclicity (checking 
whether a cycle has been already created) takes the same amount of time by 
looking for a pair of nodes in {ci,... Cj} x',y' such R{x', y') = R{y', y') = 1. 

Lemma 4. Algorithm\^ outputs in 0(\E']^^\ log |Vm|) time the first legal orien¬ 
tation of C, and each of the remaining ones with 0(log^ |Uvf|) delay. 

Proof. Before calling the LegalOrientations procedure we have to compute 
the reachability matrix from scratch and this costs 0(|A)^^| log |Vm|) time. In 
the following we will bound the delay between two outputs returned by the 
LegalOrientations procedure. Firstly, note that each call produces at least one 
solution. This is true when j = 1 since we have two possible legal orientations 


Algorithm 2: Returning all legal orientations of C 

Input-. M' {Vm ,E' m) acyclic, a cjc\eC{VcjEc) with Vc C Vm 
Output: All the legal orientations C {Vc,E c) 

Build the reachability matrix R for the nodes of Vc in M' 

Let Vc = {ci,. .. Ch}, where Ch+i = ci by definition 
Execute LegalOrientations (C'{0, 0), 1, R, 0) 

Procedure LegalOrientations {Vc,~^c), j, R, 
if j — h + Ithen 

I output C' and its set / of broken edges 
else 

U Cj, V <— Cj+1 

Ri R updated by adding the arc (u, v); 

if Ri is cyclic or has pos itiv e re ach ability test on {cj+i,... , Ch+i} then 
|_ LegalOrientations (C'{Vc,EcU{{u,v)}),j + l,Ri,I) 

R 2 -f- R updated adding the arc (v,u); 

if R 2 is cyclie or has pos itiv e re ach ability test on {cj+i,... , Ch+i} then 
|_ LegalOrientations (C'(Vc, R ^ U {(v, m)}), j ' + 1, R 2 , 1) 
if {u, t} is a ehain edge then 

if R is cyclie or has positive reachability test on {cj+i ,..., Ch+i} 

then 

|_ LegalOrientations (o', j’ + 1, R, 7 U {{m, t}}) 


of C. Before performing any call at depth j, the caller function checks whether 
this will produce at least one solution. Only calls that will produce at least one 
solution are then performed. This means that in the recursion tree, every internal 
node has at least a child and each leaf corresponds to a solution. Hence the delay 
between any two consecutive solutions is bounded by the cost of a leaf-to-root 
path and the cost of a root-to-the-next-leaf path in the recursion tree induced by 
LegalOrientations. Since the height of the recursion tree is 0(log | VmI)? be. the 
edges of C, and the cost of each recursion node is 0(log^ I^mD) we can conclude 
that the delay between any two consecutive solutions is bounded by 0(log^ I^m !)• 
As it can be seen, it is crucial that the size of C is (poly)logarithmic. □ 


Lemma 5 (Correctness). 

1. All the extended eyelie orientations of M are output. 

2. Only extended cyclic orientations of M are output. 

3. There are no duplicates. 

Proof 1. Aw extended cyclic orientation A? can be seen as the union A?" of 
A7' and C, which are two edge disjoint directed subgraphs. Our algorithm 
enumerates all the extended orientations of A7', and, for each of them, all 
legal extended orientations C : if there is a cycle in M involving only edges in 
E'j^ all the extended orientations of C are legal; otherwise just the extended 
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3. 


orientations C oi C whose arcs create a cycle in M" are legal. Hence any 
extended cyclic orientation A? is output. 

Any output solution is an extended orientation: each edge in M' and in C has 
exactly one direction or is broken. Moreover, any output solution contains 
at least a cycle: if there is not a cycle in Ad' ^ a cycle is created involving the 
edges in C. 

All the extended orientations M = M" in output differ for at least an arc in 
or an arc in Ac. Hence there are no duplicate solutions. 

□ 


As a result, we obtain delay 0 (|Am|)- 

Lemma 6. The extended cyelie orientations of M{Vm,Em) can he enumerated 
with delay 0{\Em\) and space 0{\Em\)- 

Proof. Finding extended orientationsM' of AT' can be done with 0{\E'j^^\) delay. 
Every time a new Ad' has been genera^d, we apply Algorithm [5] By Lemma H] 
we output the first cyclic orientation M of Ad with delay 0{\Em \ log \ Vm\) and 
the remaining ones with delay 0(log^ I^m|)- Hence the maximum delay between 
any two consecutive solutions is 0{\E'^ \ + \Em \ log \Vm\) = 0{\Em \ log \ Vm\) = 
0{\Em\)- The space usage is linear in all the phases: in particular in Algorithmic] 
the space is 0(log^ I^mI), because of the reachability matrix R, which is smaller 
than 0 (|Am|)- □ 

Applying ?? and ??, and considering the setup cost in Section |3| (| Vm | < |1^| 
and I Am I < lAj), we can conclude as follows. 

Theorem 1. Algorithm\^lists all cyclic orientations of G(V, E) with setup cost 
0(|Fp), and delay 0(|A|). The space usage is 0(|A|) memory cells. 


5 Absorbing the setup cost 

In this section, we show how to modify our approach to get a setup time equal 
to the delay, requiring space 0{\V\ ■ |A|). 

Theorem 2. All cyclic orientations of G(V, A) can he listed with setup cost 
0(|A|), delay 0(|A|), and space usage of 0{\V\ ■ |A|) memory cells. 

We use n = \V\ and m = |A| for brevity. Let Ai be the following algorithm 
that takes Ti = 0(mn) time to generate n solutions, each with 0(rn) delay, 
starting from any given cycle of size > logn. This cycle is found by performing 
a BFS on an arbitrary node u, and identifying the shortest cycle Cu containing 
u. Note that Cu is a log-hole as required. Now, if |Gu| < logn, we stop the setup 
and run the algorithms in the previous sections setting C = Cu. The ca^ of 
interest in this section is when \Cu\ > logn. We take a cyclic orientationGii of 
Cm and then n arbitrary orientations of the edges in G \ G„. The setup cost is 
Oim) time and we can easily output each solution in 0(m) delay. We denote 
this set of n solutions by Z\. 
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Also, let A2 be the algorithm behind Theorem[TJ with a setup cost of Oijnn) 
and 0 {m) delay (i.e. Algorithm [1]) . We denote the time taken by A2 to list the 
first n solutions, including the 0 {mn) setup cost, by T2 = 0 (jnn), and this set 
of n solutions by Z2 ■ Since Zi and Z2 can have nonempty intersection, we want 
to avoid duplicates. 

We show how to obtain an algorithm A that lists all the cyclic orientations 
without duplicates with 0 {m) setup cost and delay, using 0 {mn) space. Even 
though the delay cost of A is larger than that of Ai and A2 by a constant factor, 
the asymptotic complexity is not affected by this constant, and remains 0 (m}. 

Algorithm A executes simultaneously and independently the two algorithms 
Ai and A2. Recall that these two algorithms take T1+T2 time in total to generate 
Zi and Z2 with 0 (m) delay. However those in Z2 are produced after a setup 
cost of 0 {mn). Hence A slows down on purpose by a constant factor c, thus 
requiring c(Ti + T2) time: it has time to find the distinct solutions in Zi U Z2 
and build a dictionary Di on the solutions in Zi. (Since an orientation can be 
represented as a binary string of length m, a binary trie can be employed as 
dictionary Z?i, supporting each dictionary operation in Oim) time.) During this 
time, A outputs the n solutions from Zi with a delay of c(Ti + T2)/n = 0 {m) 
time each, while storing the rest of solutions of Z 2 \ Zi in a buffer Q. 

After c(Ti +T 2 ) time, the situation is the following: Algorithm A has output 
the n solutions in Z\ with 0 {m) setup cost and delay. These solutions are stored 
in Di, so we can check for duplicates. We have buffered at most n solutions of 
Z2 \ Z\ in Q. Now the purpose of A is to continue with algorithm A2 alone, 
with 0 {m) delay per solution, avoiding duplicates. Thus for each solution given 
by A2, algorithm A suspends A2 and waits so that each solution is output in 
c(Ti +T 2 )/n time: if the solution is not in Di, A outputs it; otherwise A extracts 
one solution from the buffer Q and outputs the latter instead. Note that if there 
are still d duplicates to handle in the future, then Q contains exactly d solutions 
from Z2\Zi (and Q is empty when A — 2 completes its execution). Thus, A never 
has to wait for a non-duplicated solution. The delay is the maximum between 
c{Ti+T2)/n and the delay of A2, hence 0 {m). The additional space is dominated 
by that of Q, namely, 0 {mn) memory cells to store up to n solutions. 

We also have an amortized cost using the lemma below, where f{x) = 0 (x) 
and s = \V\. 

Lemma 7. Listing all the extended cyelie orientations of M{Vm, Em) with de¬ 
lay 0 (/(|Am|)) and setup cost 0 (s ■ |Hm|) implies that the average cost per 
solution is 0 (/(|Am|) + |Am|). 

Proof. We perform a BFS on an arbitrary node u, and identify the shortest cycle 
Cu{Vu,Eu) that contains u. This costs 0(m) time. Note that Cu{Vu,Eu) is a 
hole (i.e. it has no chords). Note that a minimum cycle in M either is Cu or 
contains a node in Vm — Vu'- hence we perform all the BFSs from each node in 
Vm — 14 , as explained in [B] with an overall cost of 0(|Hm| • \Vm — 14|). The 
number of extended orientations of M is at least Q^J. 

setup cost is 0 {s ■ IHmI), with s < \Vm\, and the number of solutions is at least 
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2'*. The overall average cost per solution is at most 


2 ^* 


= 0[fi\EM\) + \EM\-^) 


□ 


6 Conclusions 

In this paper we considered the problem of efficiently enumerating cyclic orienta¬ 
tions of graphs. The problem is interesting from a combinatorial and algorithmic 
point of view, as the fraction of cyclic orientations over all the possible orienta¬ 
tions can be as small as 0 or very close to 1. We provided an efficient algorithm 
to enumerate the solutions with delay 0{m) and overall complexity 0{a ■ m), 
with a being the number of solutions. 
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